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j^ \ Abstract 



We discuss the signatures of non-Gaussianity in the curvaton model where the 

potential includes also a non-quadratic term. In such a case the non-linearity pa- 

Q^ I rameter /nl can become very small, and we show that non-Gaussianity is then 

^O . encoded in the non-reducible non-linearity parameter ^nl of the trispectrum, which 

can be very large. Thus the place to look for the non-Gaussianity in the curvaton 

model may be the trispectrum rather than the bispectrum. We also show that (^nl 

f^ ■ measures directly the deviation of the curvaton potential from the purely quadratic 

00 . form. While ^nl depends on the strength of the non-quadratic terms relative to the 

quadratic one, we find that for reasonable cases roughly (^nl '^ 0(— 10^) — C'(— 10^), 

which are values that may well be accessible by future observations. 



1 Introduction 

The high precision of cosmo logical observations such as WMAP [1,2] and the forthcoming 
Planck Surveyor Mission [3] will soon make it possible to probe the actual physics of the 
primordial perturbation rather than merely describing it. As is well known, the primordial 
scalar and tensor power spectra, characterized by spectral indices the and tensor-to-scalar 
ratio, can be used to test both models of inflation and other mechanisms for the generation 
of the primordial perturbation such as the curvaton [4-6] or the modulated reheating 
scenario [7,8]. However, many models can imprint similar features on the primordial power 
spectrum. Therefore the possible non-Gaussianity of the photon temperature fluctuation, 
which may provide invaluable implications on the physics of the early universe, has been 
the focus of much attention recently. 

The simplest inflation models generate an almost Gaussian fluctuation. In contrast, 
in the curvaton scenario there can arise a large non-Gaussianity [9-15]. However, in most 
studies on the curvaton so far, one simply assumes a quadratic curvaton potential. Since 
the curvaton cannot be completely non- interacting (it has to decay), it is of interest to 
consider the implications of the deviations from the exactly quadratic potential, which 
represent curvaton self-interactions. Such self-interactions would arise e.g. in curvaton 
models based on the flat directions of the minimally supersymmetric standard model 
(MSSM) [16]. Even small deviations could be important for phenomenology, as was pointed 
out in [11] where it was shown that the non-Gaussianity predicted by the curvaton model 
can be sensitive to the shape of the potential. In particular, the nonlinearity parameter 
/nl which quantifies the bispectrum of primordial fluctuation can, in contrast to the case 
of the quadratic curvaton potential, be very small in some cases. 

However, the signatures of non-Gaussianity can be probed not only with the bispectrum 
but also with the trispectrum. Thus even if /nl is very small, it does not necessarily imply 
that the fluctuation is almost Gaussian; rather, in that case the imprint of non-Gaussianity 
may be detected only in higher order statistics. The situation is then more complicated 
since trispectrum is characterized by two numbers, conventionally denoted as tnl and (^nl- 
Roughly, Tnl ~ /nl ^^^ single field models including the curvaton, whereas ^^nl measures 
the four-point correlator that is not reducible to the three-point correlator. In general, 
one can expect that when /nl is large, also (7nl is large, although this is not a necessity, 
as will be discussed below. 

The main purpose of this paper is to investigate non-Gaussianity in the curvaton 
model by considering both the bispectrum and the trispectrum assuming a potential which 
includes also a non-quadratic term. We will show that even when /nl is negligibly small, 
(^NL can be very large, indicating that the first signature of the curvaton-induced primordial 
non-Gaussianity may not come from the bispectrum but rather from the trispectrum. 

The structure of the paper is as follows. In section 2, we summarize the formalism, 
describe the set-up of the curvaton model, and give some key formulae for the study of 
non-Gaussianity in the curvaton model. Then in section [31 we investigate non-Gaussianity 
in curvaton models, presenting the predictions for the nonlinearity parameters /nl and 



5'NL- We also discuss how the predictions are affected by the size of the couphng and the 
power of a non-quadratic term. The final section is devoted to the conclusions. 

2 Non-linearity parameters 

Our starting point is the 6N formalism [17-20] that can be used to calculate the primordial 
curvature perturbation (. In this approach, the final amplitude of the primordial curvature 
perturbation on the uniform energy density hypersurface is given by 6N{t, x) where 6N 
is the perturbation of the number of e-folds on a uniform energy density hypersurface 



during the final radiation dominated epoclcj with respect to the the initial fiat time-slice 



at horizon crossing during inflation. In the following, we are interested in the fluctuations 
generated from the curvaton field; thus we can write ( as 

C = NJa, + -N,,{6a,Y + -K,,{6a,Y + ■■■ (1) 

2 o 

where a is the curvaton field and the derivative of A^ with respect to cr are represented as 
N^ = dN/da, N„„ = d'^N/da'^ and N,,, = d^N/da\ 

To discuss non-Gaussianity in the scenario, we make use of the non-linearity parameters 
/nl and ^^nl defined by the expansion 

C = Ci + ^/nlCi' + ^9nlC! + 0{Ct). (2) 
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Writing the power spectrum as 

(Cfc,CA;,) = (2^)'n(^^i)^(^i + ^2), (3) 

the bispectrum and trispectrum are given by 

{CuMU) = i^^fBdku ^2, ks)5ih + k2 + h). (4) 

(Cfc^CfeCfcaCfc,) = WT^h, h, h, h)S{h + h + h + h), (5) 

where S^ and T^ are products of the power spectra and can be written as 

Bciku k2, h) = ^/nl (Pc(^i)n(^2) + Pdk2)Pdk3) + Pc{k3)Pc{ki)) , (6) 

T^ki, /C2, h, ki) = tnl (Pc(^13)Pc(^3)Pc(^4) + 11 perms.) 

54 
+^^NL (Pc(^2)Pc(^3)Pc(^4) + 3 perms.) . (7) 



'^-^ In the curvaton scenario, there may arise two radiation dominated epochs. The final radiation epoch 
here means the one after the curvaton decay. 



Note the appearance of the independent non-hnearity parameter tnl- However, in the 
scenario we are considering in the following, tnl is related to /nl by 



Tnl 



36 
25 



/nl 



To evaluate the primordial curvature fluctuation in the curvaton model, we need to 
specify a potential for the curvaton. Here we go beyond the usual quadratic approximation 
and consider the following potential: 



V{a) 



m.<7 + \m^ 



a 



(9) 



which contains a higher polynomial term in addition to the quadratic term. For later 
discussion, we define a parameter s which represents the size of the non-quadratic term 
relative to the quadratic one: 

Thus the larger s is, the larger is the contribution from the non-quadratic term. 

When the potential is quadratic, the fluctuation evolves exactly as the homogeneous 
mode. However, when the curvaton field evolves under a non-quadratic potential, the 
fluctuation of the curvaton evolves non-linearly on large scales. In that case the curvature 
fluctuation can be written, up to the third order, as [13] 
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where cTosc is the value of the curvaton at the onset of its oscillation while r roughly 
represents the ratio of the energy density of the curvaton to the total density at the time 
of the curvaton decay. The exact definition is given by 



3/0^ 
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Notice that a'^^^/aosc = l/c* for the case of the quadratic potential. With this expression, 
we can write down the non-linearity parameter /nl as 



/nl = ^I1 



cr: 
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(13) 
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Figure 1: (Left) Plot of /nl as a function of n for several values of s. (Right) Plot of 
(7nl as a function of /nl for several values of s. Notice that /nl and n have one-to-one 
correspondence. In both panels, r = 0.01. 



Also notice that, although the curvaton scenario generally generates large non-Gaussianity 
with /nl ^ C^(l)) the non-linearity parameter /nl can be very small in the presence of the 



non-linear evolution of the curvaton field which can render the term 1 
[11,13]. 

In the curvaton model, the non-linearity parameter gfNL can be written as 
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(14) 
As one can easily see, even if the non-linear evolution of o cancels to give a very small 
/nl, such a cancellation does not necessarily occur for gfNL- This indicates an interesting 
possibility where the non-Gaussian signature of the curvaton may come from the trispec- 
trum rather than from the bispectruni^^l. In the next section, we discuss the ramifications 
of this possibility in detail. 



3 Signatures of non-Gaussianity 

Let us now consider non-Gaussianity in curvaton models, paying particular attention to 
the non-linearity parameters /nl and (^nl- We have studied their behavior numerically 
and plotted the results in Figs. dHU 

In the left panel of Fig. [H the value of /nl is plotted as a function of the power n for 
several values of s. There we have fixed the value of r to r = 0.01. As can be read off 



"^^ Some general discussions on this issue can be found in Rcf. [13]. 





Figure 2: (Left) Plot of the combination (7oscC"osc/^osc ^s a function of n for several values 
of s. (Right) Plot of the combination cr^g^a^'g^/cr^g^ as a function of n for several values of 



s. 



from Eq. (IT3|) . when r is small, /nl can be well approximated as 






cr. 



/2 



(15) 



Notice that the combination o"oscC"osc/^osc i^ ^lero when n = 2. As the value of n becomes 



larger, the above combination yields a negative contribution, which is depicted in the left 
panel of Fig. [2j Thus at some point, (o"oscO"osc)/c"osc becomes ~ —1 and gives /nl ~ 0. In 
other words, /nl decreases to zero as the potential deviates away from a quadratic form 
and then becomes zero for some values of n and s (see also [11]). To see in what cases 
we obtain /nl = 0, we show the contours for /nl in Fig. O As seen from the figure, for 
small values of s, which correspond to the cases where the size of the non-quadratic term 
is relatively small compared to the quadratic one, the power n should be large to make 
/nl very small. It should also be mentioned that, for a fixed s, if we take larger values of 
n beyond the "cancellation point" /nl = 0, /nl becomes negative. 

However, as already discussed, even if we obtain very small values for /nl, it does not 
necessarily indicate that non-Gaussianity is small in the model but may show up in the 
higher order statistics. Indeed, this appears to be a generic feature of the curvaton model: 
the trispectrum cannot be suppressed and (^nl can be quite large and is always negative 
for small values of r. This can be seen in the right panel of Fig. [H where we plot the value 
of (7nl as a function of /nl- Note that once we choose the value of n, /nl is given for fixed 
s and r. Thus a function of /nl can be regarded as a function of n. The corresponding 
value of n for /nl can be read off from the left panel of Fig. [T] (We also show the plot of 
(y^NL as a function of n in the top left panel of Fig. HI) 

Interestingly, even if /nl is zero, the value of Is'nlI can be very large, as can be seen 
Fig. [H This is because a cancellation which can occur for /nl does not take place for (yfNL, 
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Figure 3: Contours of /nl are shown for the case r = 0.01 in the ns plane. Note that the 
hne of /nl = is not changed for small values of r. 



which is a smooth function of n. When r is small, (^nl is mainly determined by the first 
term in Eq. ([1^ 



fi-NL ^ rr 
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(16) 



In the right panel of Fig. [21 we plot the combination of (o'osc^oD/'^osc ^^ a function of n for 
several values of s. There one sees that both terms in Eq. (fT6|l give negative contributions, 
which indicates that the cancellation between these terms never occurs. Therefore (^nl is 
always negative for small values of r even when /nl is very small. Thus, it may turn out 
that the best place to look for non-Gaussianity in curvaton models is the trispectrum and 
S^NL in particular. Moreover, since for the quadratic potential with small r, /nl — 5/4r and 
(7nl — — 10/3r, in the purely quadratic case there is a relation between these parameters 
given by 

9nl ^ -y/NL • (17) 

Thus any deviation from the relation Eq. fITTj) can indicate that the potential is not 
quadratic. Therefore the non-linearity parameter (yfNL can be a direct measure for the 
deviation of the curvaton potential from the quadratic one and hence probes directly the 
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Figure 4: (Top left) Plot of — (^nl as a function of the power n for several values of s. The 
value of r is fixed as r = 0.01 in the top panels. (Top right) Plot of —g^L as a function 
of s for several values of the power n. (Bottom left) Contours of (^nl in the n-r plane for 
the case s = 0.1. (Bottom right) Contours of (^nl in the n-r plane for the case s = 0.01. 



underlying physics of the model. The larger the deviation from the purely quadratic case, 
the more negative (^nl is- Note that the right panel of Fig. [T] is somewhat confusing in 
this respect. To see how g^i^ depends on the parameters, we show several plots for gfNL in 
Fig. m In the top left and right panels of Fig. H] we plot (^nl for different values of s and 
n for a fixed r = 0.01, respectively, whereas in the bottom panel of the figure, contours of 
g'NL in the n-r plane for fixed s = 0.1 (bottom left) s = 0.01 (bottom right). 



4 Conclusion 

We have discussed the signatures of non-Gaussianity in the curvaton model with the 
potential including a non-quadratic term in addition to the usual quadratic term. When 
the curvaton potential is not purely quadratic, fluctuations of the curvaton field evolve 



nonlinearly on superhorizon scales. This gives rise to predictions for the bispectrum, 
characterized by the non-hnearity parameter /nl, which can deviate considerably from 
the quadratic case. In particular, depending the power and the relative strength of a 
non-quadratic term, the value of /nl can become zero, as already was pointed out in 
Ref. [11]. In this paper, we investigated this issue by considering also the trispectrum, 
paying particular attention to the non-linearity parameter (^nl which quantifies the the 
non-reducible part of the trispectrum. The second non-linearity parameter describing the 
trispectrum, denoted as tnl, is proportional to /^l and hence is small whenever /nl is 
small. In contrast, by studying /nl and (^nl simultaneously, we find that even when /nl 
is negligibly small, the absolute value of (^nl can be very large. Thus the signature of 
non-Gaussianity in the curvaton model may come from the trispectrum rather than from 
the bispectrum if its potential deviates from a purely quadratic form, as one would expect 
in realistic particle physics models. 

Moreover, we have shown that the deviation of (^nl from the relation Eq. flT7|) . valid for 
the purely quadratic case in the small r limit, is a direct measure of the deviation from a 
quadratic curvaton potential. In the small /nl limit, typical values of s'nl depend on the 
strength of the non-quadratic terms relative to the quadratic one, denoted as s in Eq. fITOl) . 
but roughly we find that (^nl ~ C(— 10^) — (9(— 10^), as can be seen in Fig. [1] and Figs. 
m Such values may well be accessible in future experiments such as the Planck Surveyor 
Mission. Hence we may soon be in a position to test curvaton models in a meaningful way 
by exploring physics beyond the simple phenomenological quadratic potential. 
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